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ABSTRACT 



The intrinsic sizes of the field galaxies with / < 26 in the Hubble and ESO-NTT Deep Fields 
are shown as a function of their redshifts and absolute magnitudes using photometric redshifts 
derived from the multicolor catalogs and are compared with the CDM predictions. Extending to 
lower luminosities and to higher z our previous analysis performed on the NTT field alone, we 
find that the distribution of the galaxy disk sizes at different cosmic epochs is within the range 
predicted by typical CDM models. However, the observed size distribution of faint (Mb > —19) 
galaxies is skewed with respect to the CDM predictions and an excess of small-size disks (Rd < 2 
kpc) is already present at z 0.5. The excess persists up to z ^ 3 and involves brighter 
galaxies . Such an excess may be reduced if luminosity-dependent effects, like starburst activity 
in interacting galaxies, are included in the physical mechanisms governing the star formation 
history in CDM models. 



Subject headings: galaxies: fundamental parameters - galaxies: evolution - galaxies: formation 
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1. INTRODUCTION 

The understanding of galaxy formation has recently undergone an appreciable progress. Observationally, 
this is due to photometric and spectroscopic information in deep galaxy fields. A corresponding theoretical 
progress has been achieved by the developement of semi-analytical approaches including the gas cooling and 
star formation processes into the well developed hierarchical clustering theories for dark matter (DM) halos; 
this allows to relate the galaxy DM circular velocities to observable, luminous properties. In this context, 
the TuUy-Fisher (TF) relation represents a typical test for the present theoretical models, as it relates the 
total luminosity of a disk galaxy to its halo circular velocity. 

However, from an observational point of view, the measure of circular velocity is limited to bright nearby 
spirals and a few very bright galaxies at intermediate redshifts (e.g. Vogt et al. 1997). Extending the TF 
relation to fainter/distant spirals is essential to test the evolution in the L, z plane of the M/L ratio predicted 
by CDM theories. 

An alternative statistical approach to connect the luminous and dynamical properties of galaxy disks is 
based on the size-luminosity relation, once a specific model is assumed to connect size to circular velocity. 
This has the advantage of exploring the dynamical evolution over a wide range of luminosities and redshifts. 

In a previous paper (Poli et al. 1999, Paper I) we applied this novel approach to the ESO-NTT Deep 
Field (Arnouts et al. 1999) to derive morphological information for the galaxies in the field down to / = 25 
after appropriate seeing deconvolution. The derived intrinsic angular sizes were then converted into physical 
sizes adopting photometric redshifts for each galaxy in the catalog (Fontana et al. 1999a). The distribution 
of sizes in the L, z plane were compared with predictions of CDM semi-analytic models of galaxy formation 
(e.g. Cole et al. 1994, Baugh et al. 1998) complemented with the specific size-velocity relation worked out 
by Mo, Mao, & White (1998) for rotationally supported disks. This analysis showed an excess of small-size 
low luminosity galaxies at small-intermediate redshifts. However, the sample magnitude limit did not allow 
an assessment of this excess at higher redshifts. Here we want to extend the study to higher z and lower L 
using the available data in the Hubble Deep field North where morphological information on the faint galaxy 
sample is available in the literature (Abraham et al. 1996; Odewahn et al. 1996; Driver et al. 1998; Marleau 
& Simard 1999). This will allow to assess if the excess observed at intermediate redshifts is an evolutionary 
effect or if it is present also at higher z, indicating the presence of remarkable physical processes not included 
in the standard CDM models. 
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2. THE GALAXY CATALOG 

The morphological information on the galaxies in the Hubble Deep Field north was obtained using the 
galaxy catalog by Marleau & Simard (1999) where the structural parameter values were derived from the 
intensity profile fitting. Exponential profiles were assumed for the disk component and a final characteristic 
disk radius in arcsec was computed for all the galaxies up to / ~ 26. 

The HDF-N galaxy catalog has been joined with the NTT Deep Field catalog used in Paper I limited 
to / ~ 25. 

To each galaxy in the joined catalog, a photometric redshift estimate was assigned with the same best 
fitting procedure applied in Paper I and in Giallongo et al. 1998. This was obtained through a comparison of 
the observed colors with those predicted by spectral synthesis models (Bruzual & Chariot 1996) including UV 
absorption by the intergalactic medium and dust reddening. The aperture magnitudes used to estimate the 
galaxy colors and the I total magnitude for each galaxy were extracted from the catalog by Fernandez-Soto 
et al. (1999). The catalog of the photometric redshifts is presented in Fontana et al. 1999a. The resulting 
typical redshift accuracy is A2; ~ 0.06 up to 2; ~ 1.5 and A2; ~ 0.3 at larger redshifts. 

As derived in other fields of similar depth, the bulk of the galaxies is concentrated at intermediate 
redshifts 2 ~ 0.5 — 1 with a tail in the distribution up to 2; ~ 6. 

In order to veriiy the effects of the background noise on the measured sizes of the faint galaxies in the 
HDF we have performed a set of simulations specifically designed to reproduce the typical conditions of real 
data. The same test has been performed in Paper I for the galaxies in the NTT Deep Field. 

The intensity profiles were reproduced as in the observed HDF images assuming the same pixel sampling. 
Assuming an intrinsic exponential profile, a series of synthetic images were constructed using different half 
light radii ranging from rhi = 0.1" to rhi = 0.9" with a step of 0.1". An average of 25 random objects were 
computed for each radius, assuming a total magnitude of 7 ~ 26 which is the limiting magnitude of the 
morphological galaxy sample. The background subtracted image of a bright star was selected in the field 
to reproduce the instrumental PSF. Its normalized profile was then convolved with the synthetic images of 
the disk galaxies. The convolved two dimensional profiles were randomly inserted in regions of the HDF far 
from very bright objects to reproduce the observed HDF galaxies with the appropriate pixel size and noise 
levels. Finally, the multigaussian deconvolution technique was applied to the synthetic data as in Paper I. 

We first notice that there is no selection bias against galaxies with large size {rhi ~ 1 arcsec) and low 
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surface brightness down to / ~ 26 since all the synthetic objects were detected. The results are shown 
in Fig. 1, where the error bars represent the dispersion around the mean due to noise in the background 
subtraction. A good match between the intrinsic and measured half light radii was obtained up to rhi ~ 0.7 
arcsec. For larger values, a slight underestimate in the measured values appears at the sample limiting 
magnitude. In any case it can be seen that, even for the faint galaxies with I ~ 26, the overall correlation 
between intrinsic and measured half-light radii is preserved in such a way that an intrinsically large, faint 
object, e.g. with r/j; ^ 0.7", can not be detected as a small sized one, e.g. with r/^i ^ 0.1". The simulation 
shows that the fraction of small size galaxies present in the HDF morphological catalog is real and is not 
due to intrinsically larger objects which have been shrunk by noise effects. 



3. THE DISTRIBUTION OF THE GALAXY SIZES IN LUMINOSITY AND REDSHIFT: 

A COMPARISON WITH CDM MODELS 

We have computed the disk linear size Ra and the absolute blue magnitude for each galaxy in the catalog 
using the color estimated redshifts as discussed in the previous section. 

In Fig. 2 we plot the distribution of the observed sizes for the HDF galaxies as a function of luminosity 

in four different redshift intervals. The filled circles represent HDF spirals with bulge-to-total ratio in the 
range 0.05 < B/T < 0.75 while asterisks represent galaxies with B/T < 0.05, most of which of irregular 
morphology (Marleau & Simard 1998; Schade et al. 1998). HDF galaxies with B/T > 0.75 are excluded 
since they are bulge dominated systems. The NTTDF galaxies of Paper I are also shown as empty squares. 

We also show in the shaded area the prediction of our rendition of the standard semi-analytical CDM 
models. This relates the luminous properties of galaxies to their circular velocity including the hierarchical 
merging of dark matter halos, the merging of galaxies inside the halos, the gas cooling, the star formation and 
the Supernovae feedback associated with the galaxies. Finally, the circular velocity of the halos is connected 
to the disk scale length using the Mo, Mao & White (1998) model; the latter correlation depends on the 
dimensionless angular momentum A whose lognormal distribution p{X) is given in Mo et al. (1998). The 
shaded area in Fig. 2 corresponds to that allowed for 0.025 < A < 0.1, the values corresponding to the 
10% and 90% points of p(A). The solid line corresponds to A = 0.05, the 50% point of p(A). The full p{\) 
distribution is taken into account in the differential size distribution of galaxies with / < 26 (normalized to 
the total number) shown in Fig. 3 for different redshift bins. A tilted CDM power spectrum of perturbations 
with n = 0.8 in an = 1 Universe with Hg = 50 km/s/Mpc has been used. 
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The full details of our semi-analytic model are given in Appendix A of Paper I together with the adopted 
set of star formation and feedback parameters. The latter set was chosen as to optimize the matching to 
the local I-band TuUy-Fisher relation for bright galaxies and the B-band galaxy luminosity function. Note 
that, since the disk velocity is ~ 20% higher than that of the DM, a small offset (~ 0.5 mag) between the 
predicted and the observed TuUy-Fisher relation is present (see Fig. 9 of Paper I) . 

Fig. 2 shows that at ^ < 1 and for faint magnitudes {Mb > —19), the observed sizes tend to occupy 
preferentially the small size region below the 50% locus of the angular momentum distribution. Correspond- 
ingly, Fig. 3 shows the excess of small {Rd < 2 kpc) size galaxies with respect to the CDM predictions. 
Indeed, for Mb > —19, the predicted average disk size is 2.1 Kpc while the observed one is 1.4 Kpc. These 
results are similar to those presented in Paper I, although extended down to Mb < —16 and with a larger 
statistics. 

The excess becomes less evident at brighter magnitudes, in agreement with recent studies (Lilly et al. 
1998, Simard et al. 1999) which indicate little evolution in the morphological properties of bright spirals in 
the CFRS sample up to z ^ 1. At z ^ 1, the larger statistics available with the present sample (with respect 
to that used in Paper I) shows that the excess persists and involves brighter galaxies {Mb ^ —20) with an 
observed average Rd — 1.3 kpc respect to the predicted Rd — 1.9 kpc. In addition, the excess appears for 
all the galaxies in the sample regardless of their morphological classification and so does not depend on the 
selection procedure adopted for the spiral sample. Furthermore, the above excesses cannot be due to the 
offset (only 0.5 mag) between the observed and the theoretical TuUy-Fisher relation, as confirmed by the 
good agreement of the Rd — Mb relation at the bright end. 

In summary. Figs. 2 and 3 indicate that small-size galaxies appear smaller and/or brighter than predicted 
by CDM at all z (indeed, even more at high z). Within the framework of the adopted standard scenario 
of disk formation, which assumes the conservation of the specific baryonic angular momentum (Mo, Mao 
& White 1999), a viable solution consists in introducing a brightening of small-size galaxies. In particular, 
we note that shifting the predicted curves toward higher luminosities results in a better fit to the data. At 
z ^1 the best fitting shift is ~ 1 mag, while at larger z the best fitting shift is ~ 1.5 mag (Fig. 2). 
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4. CONCLUSIONS AND DISCUSSION 

The present analysis performed on a larger and deeper sample, confirms our previous findings a,t z ^1 
(Paper I), where an excess of faint {Mb > —19), small-size {Rd < 2 kpc) galaxies with respect to the CDM 
predictions was found. The results presented here show that the excess persits even at higher redshifts 
{1 < z < 3.5) and for brighter galaxies {Mb > —20). Several processes may be responsible for the above 
excess (see Paper I), like the non conservation of the gas angular momentum during the collapse in the Mo 
et al. model. Alternatively, within the Mo et al. framework for disk formation, a possible explanation can be 
sought in luminosity-dependent effects related to the physical mechanisms involved in star formation activity 
already at high z. Indeed, a shift of the shaded region (the CDM prediction) by ~ 1 mag at ^; < 1 and 1.5 
mag at > 1.5 is sufficient to reconcile the CDM predictions with the observations. 

Such shift could be due to the starbust brightening of the numerous small size galaxies subject to close 
encounters/interactions. This brightening would have the advantage of explaining at the same time the flat 
shape of the global cosmological star formation rate M* observed at ^; > 1.5 in deep surveys (Steidel et al. 
1999; Fontana et al. 1999b) which results a factor ^ 5 — 10 higher than predicted by CDM. Since the SFR 
is proportional to the UV-B luminosity, a brightening of the predicted luminosities in small size galaxies 
by ~ 1 — 2 mag is needed in both cases. The interaction rate needed to reconcile the CDM evolutionary 
scenario with the various observables can be derived by the following simple considerations. The SFR in a 
galaxy halo from a cool gas of mass Mcooi can be written as M* « fgas Mcooi/ti + (1 — fgas) Mcooi/t*, where 
fgas is the fraction of gas converted in stars due to interactions, Tj is the timescale for interactions and r* is 
the quiescent star formation time scale. While only the latter term is usually included in the CDM models, 
we note that for fgas ~ 0.1 a r, shorter than r* by about a factor 100 would be implied to obtain a SFR 
consistent with the high z data; for the population of small galaxies (dominating at high z) with a circular 
velocity Vc <, 100 km/s (corresponding to r* 5 Gyr, see Cole et al. 1994) this would imply n « 0.1 Gyr, 
in fact close to the dynamical time scale of these systems at z ^ 2. Note that for larger {vc > 200 km/s) 
disk galaxies, t^, ^ v^^'^ becomes smaller while ^ ^/Ng remains large due to their small number density 
Ng (Cavaliere & Vittorini 1999), so that the quiescent star formation mode prevails for these systems; in 
addition since Ng r-^ {1 + z)^ for all the galaxies, the interaction time scale Tj becomes ineffective at small z. 

Detailed implementation of the interaction-driven star formation mode in semi-analytical models will 
soon provide a more quantitative test for the importance of this physical mechanism in determining the 
galaxy properties at high z. 
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FIGURE CAPTIONS 



Fig. 1. Deconvolved half light radii as a function of true values in simulated data of the HDF. Error bars 
are one sigma confidence intervals. 

Fig. 2. Distribution of galaxies in the luminosity-size plane in four redshift intervals. The disk radii are in 
kpc. Empty squares are NTTDF galaxies; filled circles are HDF spiral galaxies with 0.05 < B/T < 0.75; 
stars are galaxies with B/T < 0.05 most of which with irregular morphology. The shaded region represents 
the region allowed by the model. The upper and lower lines correspond to the 90% and 10% points of the 
angular momentum distribution. The solid line corresponds to the 50% point of the same distribution (see 
the text for details). 

Fig. 3. Size distribution of the low and high luminosity spiral galaxies in the HDF and NTTDF shown in 
Fig. 2. The corresponding curves are the distributions predicted by the CDM model. An excess of small size 
galaxies with respect to the CDM predictions is apparent at i?^ ~ 1 kpc. 
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